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ABSTRACT 

We investigate departures from LTE in the line formation of Fe for a number of well-studied 
late-type stars in different evolutionary stages. A new model of Fe atom was constructed from 
the most up-to-date theoretical and experimental atomic data available so far. Non-local ther- 
modynamic equilibrium (NLTE) line formation calculations for Fe were performed using ID 
hydrostatic MARCS and MAFAGS-OS model atmospheres, as well as the spatial and tem- 
poral average stratifications from full 3D hydrodynamical simulations of stellar convection 
computed using the Stagger code. It is shown that the Fe I/Fe II ionization balance can be 
well established with the ID and mean 3D models under NLTE including calibrated inelastic 
collisions with H I calculated from the Drawin's (1969) formulae. Strong low-excitation Fe I 
lines are very sensitive to the atmospheric structure; classical ID models fail to provide consis- 
tent excitation balance, particularly so for cool metal-poor stars. A better agreement between 
Fe I lines spanning a range of excitation potentials is obtained with the mean 3D models. 
Mean NLTE metallicities determined for the standard stars using the ID and mean 3D models 
are fully consistent. Also, the NLTE spectroscopic effective temperatures and gravities from 
ionization balance agree with that determined by other methods, e.g., infrared flux method 
and parallaxes, if one of the stellar parameters is constrained independently. 
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1 INTRODUCTION 

Accurate determination of basic stellar parameters is fundamental 
for calculations of chemical composition, ages, and evolutionary 
stages of stars. One of the most commonly used methods to deter- 
mine effective temperature T e g , surface gravity log g and metallic- 
ity [Fe/H], is to exploit excitation-ionization equilibria of various 
chemical elements in stellar atmospheres. Iron, with its partly filled 
3d subshell, has, by far, the largest number of lines all over the 
spectrum of a typical late-type star. This atomic property coupled 
to a relatively large abundance makes it a reference element for 
spectroscopic estimates of stellar parameters. 

The goal of this work is to study systematic uncertainties in 
this method, which are related to 1 ) the assumption of local thermo- 
dynamic equilibrium (LTE) in spectral line formation calculations, 
and 2) the use of theoretical ID hydrostatic model atmospheres. 
These two approximations are inherent in most of the line forma- 
tion codes utilized in spectroscopic studies, since they strongly re- 
duce the complexity of the problem permitting the analysis of very 
large stellar samples in short timescales. Yet, in conditions when 



E-mail : mbergema @ mpa-garching . mpg.de 



the breakdown of LTE/ ID modelling occurs the inferred stellar pa- 
rameters may suffer from large systematic biases (e.g., Asplund 
2005). To assess the extent of the latter, more physically realistic 
modeling is necessary. 

Non-local thermodynamic equilibrium (NLTE) effects on the 
Fe i/Fe n level populati ons for FGK stars have been extensively dis- 
cussed in the literature ( Athay & Lites 1972: Cram, Lites & Ruttenl 
1 1980 1: iBikmaev. Bobritskii & Sakhibullinll 1990l: iTheyenin & Idiard 
ll999l : lGehren et alj|200ll : iGehren. Korn & Shill200ll ; Collet et al. 
2005; Mashonkina et al. 2011). These and other studies showed 
that NLTE effects in the ionization balance of Fe i/Fe n are large for 
giants and metal-poor stars. The effect on solar-metallicity stars is 
smaller, but it is must be taken into account if one aims at the accu- 
racy of few percent, as is the case for the Sun. Despite major efforts 
aimed at understanding how non-equilibrium thermodynamics af- 
fects the line formation of Fe, there have been only few attempts 
to quantify these deviations in a systematic m anner across the 
Hertzsprung-Russell diagram. lMashonkinal d201 ll) provided a small 
grid of NLTE corrections to five Fe i lin es for the solar-metallicity 
stars with T cff > 6500 K and logg > 3. iThevenin & Idiard ( fl999l) 
explored a larger range of stellar parameters including FGK stars 
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down to [Fe/H] « -3, however, they did not present a regular grid 
of NLTE corrections. 

Furthermore, it has only recently become possible to perform 
full time-dependent, 3D, hydrodynamical simulations of radiative 
and convective energy transport in stellar atmospheres. Such sim- 
ulations (Nordlund & Drawins 1990, Asplund et al. 1999, Collet 
et al. 2006, Ludwig et al. 2006. iNordlund. Stein & Asplund|[2009l 
Freytag et al. 2010) have evidenced important shortcomings of ID, 
stationary, hydrostatic models. Especially at low metallicity, it has 
been realized that ID models, by necessarily enforcing radiative 
equilibrium, overestimate the average temperatures of shallow at- 
mospheric layers, with profound implicatio ns f or the spectral line 
formation jCollet. Asplund & Trampedacbl e.g.. 1 20071) . 

In this work, we perform NLTE line formation calculations 
for Fe i and Fe n using ID hydrostatic and mean 3D model atmo- 
spheres obtained from temporal and spatial averaging of 3D hydro- 
dynamical simulations (hereafter, (3D)). A new model of Fe atom 
is constructed from the most up-to-date theoretical and experimen- 
tal atomic data available so far. Non-local thermodynamic equilib- 
rium Fe abundances, effective temperatures, and surface gravities 
are derived for the Sun, Procyon, and four metal-poor stars. The 
efficiency of thermalization caused by inelastic H i collisions is 
calibrated as to satisfy ionization equilibrium by scaling the classi- 
cal Drawin (1968, 1969) formulae. In the next paper in the series 
( iLind. Bergemann & Asplunoj2012L hereafter Paper II), we discuss 
NLTE effects on Fe i and Fe n over a wide range of stellar param- 
eters. That paper also presents a large grid of ID NLTE abundance 
corrections for a wealth of lines in metal-rich and metal-poor dwarf 
and giant spectra. 

Before proceeding with the description of the methods, we 
shall point out a few important aspects of our study. First, due to 
a comparative nature of the analysis (ID LTE vs. ID NLTE and 
(3D) NLTE) no attempt is made to fine-tune various parameters in 
order to achieve full agreement with other results in the literature. 
Second, although by the use of the averaged 3D models we roughly 
account for hydrodynamic cooling associated with convective over- 
shoot in the simulations, the effect of horizontal inhomogeneities is 
not addressed because such calculations with our new realistic ex- 
tensive model atom are beyond current computational capabilities. 
We note, however, that our results are expected to closely resem- 
ble any future full 3D NLTE calculations once these become fea- 
sible. The reason is that in NLTE the Fe i line formation is largely 
dictated by the radiation field (as explained in detail below) orig- 
inating in deep atmospheric layers, where the significance of the 
atmospheric inhomogeneities is greatly reduced. Detailed tests of 
the (3D) models, which involve comparison with other observable 
quantities, will be presented elsewhere (Collet et al., in prep.). Full 
3D NLTE calculations with a simpler model atom of Fe will be a 
subject of a forthcoming publication. 

The paper is structured as follows. Section[2]gives a brief de- 
scription of the input model atmospheres, model atom, and pro- 
grams used to compute NLTE populations and line formation. The 
results of statistical equilibrium calculations are presented in Sec.[3] 
The analysis of the solar spectrum along with the re-determination 
of the Fe abundances for the Sun and Procyon is described in detail 
in Sec. [4] Section POl presents and discusses metallicities, temper- 
atures, and gravities obtained for the metal-poor stars. Finally, a 
comparison with stellar evolution calculations is given in Sec. 14.61 
A short summary of the work and conclusions are given in Sec. [5] 



2 METHODS 



2.1 Model atmospheres 



In the present stud y, we employ both mafag s-os dGruppl 
l2004al iGruppI l2004bl iGrupp. Kurucz & Tanl |2009|) and mar cs 



l Gustafsson. Edvardsson. Eriksson. Jergensen. Nordlund & Pled 
120081) ID LTE hydrostatic model atmospheres, as well as 3D hy- 
drodynamical models computed with the Stagger code (Galsgaard 
& Nordlund 1995, Collet et al. 201 1 ). 

The mafags-os models for all investigated stars were kindly 
provided by F. Grupp. The models are plane-parallel with convec- 
tive energy transport based on the turbulent convection model of 
ICanuto & Mazzitellil ( fT99lh . Compressible turbulence is accounted 
for using a mixing length, a cm = l/nfi\= 0. 82. The reasons for 
this choice are discussed in iBernkopj 1 1998[). We note that this 
a cm is also consistent with the results of lFrevtag. Ludwig & Steffer] 
dl999h . who attempted to calibrate the mixing length parameter us- 
ing the 2D radiation hydrodynamics simulations of stellar convec- 
tion. Line absorption is computed by the method of opacity sam- 
pling including ~ 86000 wavelength points from extreme UV to 
far IR. Extensive line lists were extracted from Kurucz and VALD 
databases. The reference solar abundances were compiled from var- 
ious literature sources, givi ng preference to NL TE determinations 
by the Munich group (e.g., iGehren et alj|200lh . In particular, the 
solar Fe abundance is set to logs = 7.50. 

The marcs models are based on opacity sampling with ~ 
10 4 — 10 5 wavelength poin ts using Kurucz and VALD linelists. Con- 
vection is included in the lHenvev. Vardva & Bodenheimej d 19651) 
version of the mixing length theory with the mixing length param- 
eter set to (2mlt = 1.5. A detailed description of the models can be 
found on the MARCS web-sitt0. For the majority of the elements, 
the reference solar abundances are that of Asplund et al. (2005). 

Three-dimensional radiation-hydrodynamic simulations of 
convection at the surface of the Sun and re ference stars were com- 
puted using the Copenhagen Stagger code dNordlund & Galsgaard 
1 1995b . The physical domains of the simulations cover a repre- 
sentative portion of the stellar surface. They include the whole 
photosphere as well as the upper part of the convective lay- 
ers, typically encompassing 12 to 15 pressure scale heights 
vertically. Horizontally, they extend over an area sufficiently 
large to host about ten granules at the surface. The simulations 
use realistic input physics, including state-of-the-art equation-of- 
state, opacities, and treatment of non-grey radiative transfer. The 
adopted refere nce solar chemical composition for the simula- 
tions is that of lAsplund, Grevesse. Sauval & Scottl j2009h . For a 
more detailed description of th e simulations' setup, we refer to 
ICollet, Magic & Asplund d201 ll) . 

For the purpose of the current study, we computed spatial and 
temporal averages of these simulations over surfaces of equal op- 
tical depth log T5000 at the chosen reference wavelength at 5000 A. 
The independent thermodynamic variables, gas density p and inter- 
nal energy per unit mass e, as well as the temperature T are first 
interpolated with cubic splines for all columns in the full 3D struc- 
ture to the reference optical depth scale. The reference scale was 
constructed to cover the relevant range for line-formation calcu- 
lations (-5< logT 5 ooo £j2) evenly in logT 50 oo- For density and in- 
ternal energy, a logarithmic interpolation is used. Other physical 
quantities, namely gas pressure P gas and electron number density 



Hp is the pressure scale height and / the so-called 'mixing length' 
http://marcs.astro.uu.se/index.php 
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Figure 1. Comparison of the temperature and gas pressure structures in the marcs and (3D) model atmospheres. Arrows indicate the optical depths, above 
which the marcs models were extrapolated. 



n e i, are looked up from the simulations' equation-of-state tables as 
a function of density and internal energy. Finally, mean (3D) strat- 
ifications are constructed by averaging the various physical quan- 
tities pertinent to line-formation calculations, lnp, T, In P gis , and 
lnn e i, on surfaces of equal optical depth logTsooo and over time 



(i.e., over all simulation snapshots). We emphasize that no hydro- 
static equilibrium was enforced after the averaging process. In the 
present application of mean (3D) stratifications to line-formation 
calculations, averages of velocity fields were not considered, and 
line broadening associated with bulk gas flows and turbulence was 
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accounted for by means of a classical, depth-independent micro- 
turbulence parameter (Sect.|2j as in ID models. 

The temperature and pressure stratifications from the ID 
marcs model atmospheres and the mean (3D) hydrodynami- 
cal models of the Sun and three metal-poor standard stars are 
shown in Fig. Q~] as a function of continuum optical depth at 
5000 A.The input parameters are the ones listed in Table Q] 
The mafags-os models are not included in the plots, because 
they adopt slightly different stellar parameters. Due to the sim- 
ilar equation-of-state, marcs and mafags-os model atmospheres 
are essentially identical in the outer layers. The differences be- 
tween ID and (3D) models are more pronounced, especially in 
the case of metal-poor stellar atmospheres. At low metallicities, 
3D stellar surface convec tion simulations predict cooler upper pho- 
tospheric stratifications (Asplund, Nordlund. Trampedach & Steinl 



1 19991 ICollet. Asplund & Trampedach! 2007 ), than corresponding 
classical, hydrostatic, stationary ID model atmospheres generated 
for the same stellar parameters. The temperature in the outer layers 
of time-dependent, 3D, hydrodynamical simulations is mainly reg- 
ulated by two mechanisms: radiative heating due to reabsorption of 
continuum radiation by spectral lines and adiabatic cooling associ- 
ated with expanding gas above granules. At low metallicities, the 
coupling between radiation and matter is weakened with respect to 
the solar-metallicity case because of the decreased line opacities; 
the adiabatic cooling term therefore prevails, causing the thermal 
balance to shift toward on average lower temperatures. Stationary, 
ID, hydrostatic models do not account for this cooling term associ- 
ated with diverging gas flows, and the thermal balance in the upper 
photosphere is controlled by radiative heating and cooling only, ul- 
timately leading, at low metallicities, to higher equilibrium temper- 
atures than predicted by 3D models (Fig.Q] left panels). The lower 
average temperatures in the outer layers of metal-poor 3D models 
also imply smaller values of the pressure scale height and, conse- 
quently, steeper pressure stratifications on a geometrical scale with 
respect to corresponding ID models. However, the lower temper- 
atures also result in lower opacities; therefore, on an optical depth 
scale, the average gas pressure in the optically thin layers of the 
metal-poor 3D models appears typically higher at a given optical 
depth than in ID models (Fig.Q] right panels). 

2.2 Statistical equilibrium codes 

The NLTE level populations of Fe i and Fe n were computed with 
an up dated version of the de tail code jButler & Giddingsl 19851) and 
MULTI dCarlssonll98o1 , ll992h . In both codes, the solution of the cou- 
pled statistical equilibrium and radiative transfer equations is ob- 
tained using an approximate lambda iteration method. In detail, the 
latter is impl emented following a f-operator ap proach of the kind 
described bv lRvbicki & Hummed dl99 lL Il992n , which allows for 
self-consistent treatment of overlapping tra nsitions and continua. 
multi is based on the method des cribed by Scharmer & C arlssonl 
dl985l) with the local operator by lOlson. Auer & Buchled Il986h . 
Scattering in the bound-bound (hereafter, b-b) transitions included 
in the model atom follows complete frequency redistribution. Up- 
per boundary conditions differ in that detail assumes no incoming 
radiation at the top, whereas multi estimates the minor contribution 
from the optically thin gas and uses a second order Taylor ex pan- 
sion of the Feaut rier variables. Ng convergence acceleration dNd 
ll974l;lAuedll987l) is implemented in both codes. 

In both codes, the background line opacity is computed with a 
Planckian source function. However, there is one difference. In de- 
tail, line opacity is consistently added at the frequencies of all b-b 



and b-f transitions in the NLTE atom. In multi, metal line opacities 
were added to the continuous opacities for the calcu l ation of pho- 
toionization rates only, as described bv lColletetai] ( l2005l) , while 
blends were neglected in the calculation of the bound-bound ra- 
diative rates. This approximation is well justified and saves com- 
putational time. Firstly, Fe itself dominates line opacity in the UV. 
Secondly, according to our tests, the upper limit to the differences 
in the NLTE equivalent widths computed with and without blends 
is ~ 0.2 percent for the Sun and 5 percent for HD 122563. 

When solving for statistical equilibrium, Fe line profiles were 
computed with a Gaussian in detail and with Voigt profiles in multi. 
Although the latter is also possible with detail, it is un-necessary. 
Our tests showed that NLTE effects in Fe i are insensitive to the 
adopted profile functions. Once the level populations were con- 
verg ed with deta il, the synthetic line profiles were computed with 
SIU dReetzl 19991) . 

Finally, we remark on the handling of the equation-of-state in 
the NLTE codes. Whereas multi includes a subroutine to compute 
ionization fractions and molecular equilibria internally, detail and 
siu require partial pressures of all atoms and important molecules 
to be supplied with a model atmosphere. These are included in 
the mafags-os models. For the marcs and (3D) models, we com- 
puted the atomic partial pressures using the equation-of-state pack- 
age built in multi. 



2.3 Model atom of Fe 

2.3.1 Levels and radiative transitions 

The number of levels and discrete radiative transitions in Fe is enor- 
mous. Recent calculations of Kuruc^] predict ~ 37 000 theoret- 
ical energy levels below and above the first ionization threshold 
of Fe i, as well as 6 025 000 radiatively-permitted transitions be- 
tween them. Such atomic models are not tractable with our NLTE 
codes; thus, we combined atomic levels and transitions into super- 
levels and super-lines (see below). In our final model, which also 
includes all experimen tal data from the NIST atomic database 
( lRalchenkoetalj2012h . the number of energy levels is 296 for Fe i 
and 112 for Fe n, with uppermost excited levels located at 0.03 
eV and 2.72 eV below the respective ionization limits, 7.9 eV and 
16.19 eV. The model is closed by the Fe m ground state. The to- 
tal number of radiatively-allowed transitions is 16 207 (13 888 Fe i 
and 2 316 Fe n). Fine structure was neglected for all levels, but the 
ground states of Fe i a 5 D (configuration ls 2 2s 2 2p 6 3s 2 3p 6 3d 6 4s 2 ) 
and Fe n a 6 D. Excitation energy of each LS term is a weighted 
mean of statistical weights and excitation energies of fine structure 
levels. The mean wavelength for a multiplet is computed from the 
weighted energy levels. 

All predicted energy levels of Fe i with the same parity above 
E low > 5.1 eV were grouped into superlevels. For the upper levels 
above 54000 cnTl (6.7 eV), we combined all levels within 1000 
cm"' (0.12 eVfl whereas below this energy limit only levels devi- 
ating by less than 10 cm"' (0.001 eV) were combined. Thus, not 
only predicted, but also some experimental levels above 5.1 eV 
were combined to superlevels. None of the Fe lines selected for 
the subsequent abundance analysis has either the lower or the up- 
per level combined. 

Transitions between the components of the superlevels 

3 http://kumcz.hai-vard.edu/atoms.html 

4 For a typical F-type star (T s g = 6000 K), the thermal energy is kT ~ 0.4 
eV. 
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Figure 2. Grotrian diagram of the Fe n (top) and Fe i atoms (bottom). Predicted levels are shown in red. Only lines selected for the abundances analysis are 
shown. 



were also grouped preserving the parity and angular momen- 
tum conservation rules. The total transition probability of a su- 
perline is a weighted average of logg/'s of individual transi- 
tions, and is computed in a nalogy to a g/-value for a multiplet 
( iMartin. Fuhr & Wieselll988l) : 

/mean = = — ^ gl X * /to (1) 

^mean 2ul §1 / u 

where / and u are indices of lower and upper levels of the un- 
grouped transition, "i mcan is the mean Ritz wavelength of a super- 
line. Grotrian diagram of the Fe model atom is shown in Fig. [2] 
Although the accuracy of individual g/-values for the multitude 



of theoretically-computed transitions is hard to quantify, we expect 
our model to be a good representation of the global atomic prop- 
erties of the Fe i/Fe n and provide physically valid description of 
statistical equilibrium of the atom in the atmospheres of late-type 
stars. 

Accurate radiative bound-free cross-sections for Fe i com- 
puted using the close-coupling method were kind ly provi ded b y 
M. Bautista (private communication 2011, lBautistal see also ll997h . 
These data are computed on a more accurate energy mesh and pro- 
vide better resolution of photoionization resonances compared to 
the older data, e.g. provided in the TOPbase. Thus, 136 levels of 
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Fe i, including different multiplet systems, from singlets to septets, 
are represented by quantum-mechanical data. Hydrogenic approxi- 
mation was used for the other levels. 



2.3.2 Collisional transitions 

The rates of transitions induced by inelastic collisions with free 
electrons (e~) and H i atoms were computed using different 
recipes. For states coupled by allowed b-b and b -f tran sitions, 
we us ed the formulas of Ivan Regemorte 3 d 19621) and ISeatonl 
dl962l). respectively, in the case of e collisions and that of 
ISteenbock & Holwegerl dl984l) . in the case of H i collisions. Also, 
all states are coupled by f orbidden tran sitions induced by e~ colli- 
sions using th e formula of lAllenl dl973l) and by H i collisions with 
the formula of lTakedal dl994h . 

We constrain the efficiency of collisions with H i empirically 
(Sect. l4l |4.51 . We apply a scaling factor Sh to the cross-sections, 
which is iteratively derived from the constraint of Fe i/Fe n ion- 
ization equilibrium for the standard st ars with independently -fixed 
stellar parameters. As pointed out by iBarklem et al.l d201lh . such 
approach is an over-simplification of the problem. However, at 
present there are no useful alternatives to this classical recipe. Ev- 
idence for the necessity of including inelastic collisions with H i 
in statistical equilibrium calculations for Fe and other elements ha s 
been demonstrated in many studies (e.g., fMashonkina et al.l201 lb . 
which is also confirmed by us (Sect.[4l l4~5t . 



3 STATISTICAL EQUILIBRIUM OF FE 

In the following section, we present the results of statistical equi- 
librium calculations for Fe obtained with detail and multi. We 
briefly describe the physical processes responsible for deviations 
from LTE level populations in the Fe i and Fe n atoms for atmo- 
spheric conditions typical of FGK stars. A comprehensive descrip- 
tion of NLTE effects as a function of stellar param eters is deferred 
to the second paper of this series dLind et al.ll2012h . The main aim 
here is to understand the differences in terms of NLTE Fe i line for- 
mation with ID and (3D) model atmospheres, which is necessary 
in order to explain the large differences between the spectroscopic 
values of T e $, log g, and [Fe/H] obtained for the reference stars in 
the two cases (Sect, |43t , 

Fe i is a minority ion in the atmospheres of late-type stars, 
which increases its sensitivity to NLTE over-ionization. The over- 
all effect is that, compared to LTE, the statistical equilibrium of Fe 
favors lower number densities of the neutral ion, Fe i, although the 
number densities of relevant Fe n levels remain nearly thermalized. 
In general, this leads to a weakening of Fe i lines compared to LTE 
that, in turn, requires larger Fe abundance to fit a given observed 
spectral line. The actual magnitude of departures and NLTE abun- 
dance corrections depends on stellar parameters (Sect. l4~5l l. 

Figure [3] presents the departure coefficients[f| of selected Fe i 
and Fe n levels computed for the marcs and (3D) model atmo- 
spheres with multi. The results are shown for the Sun, the metal- 
poor subdwarf HD 84937, the metal-poor sub-giant HD 140283 
and the metal-poor giant HD 122563. Only a few, selected levels, 
typical for the dependence of their associated departure coefficients 
with depth, are included in the plot: a 5 D 4 (ground state of Fe i), 
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Figure 4. Net bound-free radiative brackets computed for the solar marcs 
model atmosphere (top and middle panel). Optical depth is indicated in each 
panel. Red and blue dots indicate over-ionization and over-recombination, 
respectively. 



z 7 D° (2.4 eV), e 5 D (5.4 eV), 1 3 H D , and two levels, which belong 
to the Fe n ion: a 6 D 9/2 (ground state) and z 4 D° (5.5 eV). In the 
optically thin atmospheric layers these and majority of other Fe 1 
levels are underpopulated compared to LTE, £>, < 1. For all stars, 
but HD 122563, the Fe 11 number densities remain close to LTE 
values throughout the atmosphere, b, ■ x 1, and a minor overpop- 
ulation of the Fe 11 ground state develops only close to the outer 
atmospheric boundary. 

Deviations from LTE in the distribution of atomic level popu- 
lations arise because the mean radiation field, J v , at different depths 
and frequencies is not equal to the Planck function, B v |T c (t)]. For 
Fe 1, excess of the mean intensity over Planck function in the UV 
continua leads to over-ionization, which sets in at log T5000 ~ 0.2, 
i.e., as soon as the optical depth in the photoionization continuum 
of low-excitation Fe 1 levels with E ss 2 eV falls below unity (Fig. 

St). 

In the layers with logTsooo < — 1, the dominant mechanism 
is over-ionization from the Fe 1 levels with excitation energies at 
2-5 eV. The lower-lying levels, including the ground state of Fe 1, 
maintain underpopulation due to radiative and collisional coupling 
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;ure 3. Departure coefficients for the four reference stars computed with the ID marcs (left panel) and (3D) (right panel) model atmospheres. 
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Figure 5. Total number densities of Fe i/Fe n for the Sun, HD 140283, HD 84937, and HD 122563. 



with the former. Excitation balance of Fe i is also mildly affected 
by radiative imbalances in line transitions. These include radiative 
pumping by the non-local UV radiative field, as well as photon suc- 
tion driven by photon losses in large-probability Fe i transitions 
between highly-excited levels. These processes leave a character- 
istic imprint on the behavior of ^-factors in the outer atmospheric 
layers, logT 50 oo < -2. In the infrared continuum, J v < B y leads to 
over-recombination, which is very efficient for our atomic model 
with only 0.03 eV energy gap of the upper Fe i levels from the Fe n 
ground state. This is illustrated in Fig. [4] which shows net radia- 
tive bracket^] for all Fe i levels in the model atom at the depths 
logTsooo = 0.12, -1.2, and -3.86. All Fe i levels with excitation 
energy between 2 and 6 eV experience net over-ionization, and the 
loop is closed by net over-recombination to the upper levels. Note 
that only radiative rates are plotted. We also compared absolute ra- 
diative rates from multi and detail and found that they agree very 
well with each other. 

A comparison of the departure coefficients computed with ID 
and (3D) model atmospheres reveals important differences (Fig. 
[3}, which depend on stellar parameters, although their qualitative 
behavior is the same. Generally, metal-poor (3D) models show 
more thermalization close to continuum-forming layers, while in 
the outer layers Fe i atoms experience a larger degree of over- 
ionization. This can be primarily understood based on the tempera- 



6 Net radiative brackets p are defined as pR c jn c = R c in c - R; c «i, where R c j 
are rate coefficients for radiative transitions from the conti nuum c to a bound 
level i of an atom and n c are atomic level populations ( Mihalas & Athay 
1 19731) . 
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Figure 6. Mean radiation field compared to the local Planck function in the 
UV continuum computed with the ID and (3D) model atmospheres of the 
metal-poor sub-giant HD 140283. 

ture gradients, which in the (3D) stratifications are shallower in the 
inner layers but steeper in the outer layers compared to ID mod- 
els. We briefly consider the case of the metal-poor subgiant HD 
140283, for which the ID and (3D) model atmospheres are com- 
pared in Fig. [T] The (3D) model leads to smaller deviations from 
LTE at -2 < logr 5 ooo S 0. That is, for the majority of the under- 
populated Fe i levels, £>,(1D) < Z>,«3D)). Inspection of the respec- 
tive T(t) relations in Fig. [T] shows that the (3D) model is cooler 
than the marcs model at < logT 50 oo S= 1> where the UV continua 
form, but it is ~ 100 K hotter higher up. Figure [6] illustrates mean 
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angle-averaged intensities at 2200 A, i.e., at the wavelength sam- 
pled by the ionization edges of the important Fe i levels. These are 
shown for the ID and (3D) models, and are compared with the local 
Planck functions B v at each optical depth. At -1.5 < log r 500 o S 0, 
the J v - By imbalance is smaller in the (3D) model compared to 
the ID model, and, thus, over-ionization is less efficient. In this 
range of optical depths, densities in the (3D) stratification are also 
slightly larger than in the ID model (Fig.[TJ. At logTsooo S -2, the 
situation is reverse. The (3D) model is nearly 1000 K cooler and 
less dense than the ID model. The J v - B Y imbalance becomes very 
large causing significant underpopulation of the lower Fe i levels, 
while uppermost Fe i levels and the Fe n ground state and excited 
levels develop appreciable overpopulation compared to their LTE 
occupation numbers. In contrast, ID models predict nearly ther- 
malized Fe n level populations, over the whole optical depth scale. 

Figure[5]shows the ratios of Fe i/Fe n number densities for the 
four cases: marcs and (3D) model atmospheres, NLTE and LTE. 
In the LTE case, about ~ 5 — 10 percent of the element is in the 
form of Fe i above logT 50 oo ~ 0. In the (3D) models, due to their 
steeper T(t) relation, LTE recombination of Fe n to Fe i is even 
more efficient, raising Fe i/Fe n to 20 - 50 percent. Under NLTE, 
this effect is suppressed since ionization balance is primarily set by 
radiative ionization and recombination. For the metal-poor stars, 
Fe i occupies less then 1 percent of the total element abundance 
at logr 5 ooo < -2 notwithstanding much cooler outer layers of the 
(3D) models. We note, however, that the huge differences in the 
ionization balance of Fe i/Fe n in the outer layers of the models are 
not important for our NLTE abundance determinations. Most of the 
Fe i lines observed in spectra of the selected metal-poor stars form 
close to the continuum forming regions, -2 < log t 30 oo < 0. 

Extensive tests demonstrated that the behavior of departure 
coefficients is robust against the differences in terms of model at- 
mospheres and NLTE codes used for the statistical equilibrium cal- 
culations. Perhaps, the only systematic effect is that multi predicts 
slightly stronger NLTE effects than detail mainly due to the dif- 
ferences in the background opacity. Also, the comparison of the 
LTE and NLTE equivalent widths obtained with multi and de- 
tail/siu show agreement to within 1-5 percent for the lines with 
W A > 10 mA; NLTE W A 's are somewhat divergent for the weakest 
lines, with the relative (multi - detail) differences being about 10 
percent. These differences shall be taken as a formal error, which 
arises as a consequence of different implementations of the same 
basic physics in the various codes. Still, this is not of a big concern 
in our work because the calibration of the models and grid calcula- 
tions are consistently performed with same code (multi). Further- 
more, the magnitude of the NLTE effects is initially effectively con- 
trolled by a free parameter, the scaling factor S H for the efficiency 
of the H i collisions. 



4 SPECTROSCOPIC STELLAR PARAMETERS 
4.1 Method 

Spectroscopic stellar parameters, including metallicity, effective 
temperature, gravity, and microturbulence, were determined for 
each model atmosphere (marcs, mafags, (3D)) in LTE and NLTE 
using the following approach. 

In a first step, we performed full spectrum synthesis for the 
reference stars with the siu line formation code using the ID 
mafags-os model atmospheres computed with parameters described 
in Sect. 14.21 The mafags-os models were used at this stage because 



they include partial pressures for molecules, which are important 
contributors to the background opacity of cooler stars. The profiles 
of all diagnostic Fe lines computed in LTE and NLTE were visu- 
ally fitted to the observed spectra. The lines were computed with 
depth-dependent Voigt profiles taking into account various exter- 
nal broadening mechanisms. The equivalent widths of the Fe lines 
were then obtained from the best-fitting NLTE profiles, excluding 
contribution of blends. The W/s are given in Tables0and|2]of the 
Appendix and are accurate to within 1-2 percent. 

In a second step, the equivalent widths W/s were applied to 
determine spectroscopic parameters using the grids of LTE and 
NLTE equivalent widths computed with multi for marcs and (3D) 
model atmospheres. For each star, we constructed a local model 
atmosphere grid with two effective temperature points (the IRFM 
value and a value 200 K lower), two values of microturbulence (1 
and 2 km s ), and a range of plausible metallicities (A[Fe/H] up to 
±1.5 dex) with a step-size of 0.25 dex. Equivalent widths were com- 
puted for the grid models in LTE as well as in NLTE for two values 
of the H i collision efficiency parameter, Sh = 0.1 and Sh = 1. 
In order to save computational time in the construction of the local 
model grids in the (3D) case, the (3D) models were simply adjusted 
by multiplying the temperature and electron pressure at a given op- 
tical depth with the ratios of those quantities obtained from corre- 
sponding ID models. Based on our experience with scaling ID and 
3D model stellar atmospheres to different stellar parameters, this 
procedure appears accurate to first-order level. The grid was thus 
constructed for marcs and (3D) models. 

mafags-os model parameters were adjusted by assuming the 
equivalent widths to have the same sensitivity to stellar parameters 
as for marcs model, which is reasonable considering the similarity 
between the two codes. 

An iterative procedure was then applied to determine spectro- 
scopic parameters from Fe i and Fe i lines. First of all the microtur- 
bulent velocity was determined by flattening the slope of Fe n line 
abundances with reduced equivalent width, thereby circumventing 
differential NLTE effects with line strength of Fe i lines. Only for 
the most metal-poor star G 64-12, the microturbulence is poorly 
defined due to the lack of strong Fe n lines. We then followed two 
different approaches, keeping either log g or r eff fixed, and optimiz- 
ing the other parameter to establish ionization balance between the 
Fe i and Fe n lines. In order to avoid biased results, we have not 
attempted to adjust the oscillator strengths of the Fe lines, giving 
preference to the experimental data from the laboratory measure- 
ments, or their weighted means if few measurements were avail- 
able (Sect. |4~3l l. Hence, all abundance results are absolute values. 
Furthermore, we did not perform a differential stellar abundance 
analysis with respect to the Sun. First, that would conflict with the 
use of the reference stars as a calibration sample. Second, a dif- 
ferential approach would introduce yet an another source of error 
related to the usual problems of fitting the strong Fe lines in the so- 
lar spectrum, in particular the ambiguity between the effect of van 
der Waals damping and abundance on a line profile. 

The standard technique to infer effective temperatures by flat- 
tening the slope of Fe i line abundances with excitation potential 
is particularly sensitive to the choice of transition probabilities. We 
found, however, that it is only weakly sensitive to Sh, which is the 
main parameter we seek to constrain. Initial attempts to vary all pa- 
rameters simultaneously, i.e. microturbulence from line strengths, 
effective temperature by excitation balance and surface gravity by 
ionization balance, did not give a conclusive answer in terms of 
the best choice of Sh- This apparent weakness is likely related to 
the problem of multidimensionality in the parameter space, and, as 



10 M. Bergemann 



Table 1. Input parameters for the reference stars. Parallaxes n and their 
uncertainties are also given. See text. 



Star 


7*eff 


(T 


Ref. 


logs 


cr 


71 




[Fe/H] 


Sun 


5777 






4.44 










Procyon 


6543 


84 


a 


3.98 


0.02 


284.56 


1.26 


-0.03 


HD 84937 


6408 


66 


b 


4.13 


0.09 


3.83 


0.78 


-2.16 


HD 140283 


5777 


55 


b 


3.70 


0.08 


17.16 


0.68 


-2.38 


HD 122563 


4665 


80 


c 


1.64 


0.16 


4.22 


0.35 


-2.51 


G 64-12 


6464 


61 


b 


4.3* 




0.57 


2.83 


-3.12 



References: a - lAufdenberg et all Hw^ : b - ICasagrande et alj <2010t) : c - 
Casagrande (private communication) 
* log g is derived from the upper limit on the parallax 

such, shall be inherent to the method itself. In addition, we do not 
exclude remaining systematic uncertainties in the atomic and atmo- 
spheric models. It is possible that a differential approach between 
similar stars would be more successful in establishing temperatures 
based on excitation balance. We therefore reduced the dimension- 
ality of the problem by not enforcing strict excitation equilibrium 
and focusing only on the ionization balance, as described above. 
As the results will show, an adequately flat trend with excitation 
potential is anyway naturally achieved for the majority of stars, HD 
122563 being a notable exception. 

The analysis described above was applied to each star in the 
sample using the marcs, mafags, and (3D) model atmospheres, thus 
yielding four desired quantities: T s g, logg, [Fe/H], and f t . The re- 
sults are discussed in Sect. 14. 51 Table[5]gives the mean metallicities 
determined using the reference T cff and log g, and abundances av- 
eraged over the measured Fe i and Fe n lines together with their 
standard deviations are presented in Table[7] 



which includes the uncertainty on the zero point of the T c g scale, 
reddening and photometric errors. 

The IRFM value for the metal-poor giant HD 122563, T eff = 
4665 + 80 K, was kindly provided by L. Casagrande (private com- 
munication). This estimate incorporates a correction due to redden- 
ing, E{B - V) = 0.005, which was determined from the interstellar 
Na D lines detected in the UVES spectrum. Our equivalent widths 
for the 5889.95 and 5895.92 A Na i lines are 22.3 and 12.4 mA, 
respectively. 

For Pro cyon, an interferometric estimate of an gular diameter 
is available dAufdenberg. Ludwig & Kervellall2005h , which gives 
r eff = 6545 ± 83 K and combined with a very accurate parallax 
7r = 284.56 ± 1.26 milli-arcsec, surface gravity can be estimated, 
logg = 3.99 ±0.02. 

For the other stars, we determined surf ace gravities from 
the Hipparcos parallaxes dvan Leeuwenl l2007j) using the m asses 



the Hipparcos parallaxes (van Leeuwen zUU / ) using tne m asses 
estimated fr o m the trac ks of VandenBerg et all d2000h by 
IGehren et al] d2004 |2006|) and apparent bolometric magnitudes 
from ICasagrande et alj d20ip|) . The uncertainties were computed 
by mapping the uncertainty in masj^ (0.2 M Q for HD 122563 and 
G 64-12, and 0.1 M for HD 84937 and HD 140283), temperature 
and parallax into the range of possible gravities (A logg < 0.16 
dex). The only exception is G 64-12, for which the parallax is too 
uncertain, n = 0.57 ± 2.83 milli-arcsec. From comparison with 
metal-poor evolutionary tracks for reasonable masses (Fig. lilt , 
one may also conclude that within the errors of the given T s g (6464 
K) surface gravity of G 64-12 is in the range logg = 3. 8... 4. 6 
dex. We, thus adopt the value derived from the upper limit on 
the parallax, logg = 4.3 and assign a nominal erro r of 0. 3 
dex. In comparison, Nissen. Chen. Asplund & Pettinil ( 2004b ; 



iNissen. Akerman. Asplund. Fabbian. Kerber. Kaufl & Pettinil 
1 20071) and iFabbian. Nissen. Asplund, Pettini & Akermanl d20( 
estimate logg = 4.26 



4.2 Observations and reference stellar parameters 

Our reference sample consists of six late-type stars (Table [TJ, and 
includes two solar-metallicity stars (the Sun and Procyon), two 
metal-poor dwarfs (HD 84937 and G 64-12), a metal-poor subgiant 
(HD 140283), and a very bright metal-poor giant (HD 122563). 
The following observed spectra were adopted here. For the Sun, 
we used the KPNO flux spectrum (Kurucz et al. 1984). The UVES 
observations of Procyon (HD 61421), HD 84937, H P 140283, and 
HP 1 22563 were taken from the UVES-POP survey dBagnulo et"al] 
l2003h . These spectra have a slit-determined resolution of A/AA ~ 
80000 and a signal-to-noise ratio S/N ~ 300 near 5000 A. The 
UVES spectrum for G 64-12 was taken from ESO/ST-ECF Science 
Archive Facility (67.O-0554(A)). The number of Fe lines suitable 
for the spectrum synthesis at this spectral quality, is 40 (Procyon) 
to 10 (G 64-12). For comparison, we also used the Keck/HIRES 
spectra of G 64-12 with A/AA ~ 100000 and S/N ~ 500 kindly 
provided by J. Melendez (private communication). 

A crucial step in our analysis is the choice of the reference T c fj 
and log g values for the selected stars, which are used as a bench- 
mark for testing the new NLTE model atom and (3D) model atmo- 
spheres. These parameters were taken from the literature, giving 
preference to the least model-dependent methods, such as interfer- 
ometry, infra-red flux method (IRFM), and parallaxes. The adopted 
values and their uncertainties are given in the Table [T] A brief de- 
scription of these data is given below. For the fou r metal-poor stars, 
we ad opted the IRFM effective temperatures by ICasagrande et all 
d2010t) . The mean internal uncertainty of the data is about 70 K, 



4.3 Line selection 

The Fe lines for the abundance calculations were selected on a star 
by star basis, i.e., carefully inspecting the observed stellar spectra. 
We rejected lines affected by blends, strong damping wings or are 
located in the spectral windows where continuum placement is un- 
certain. Our solar line list includes 59 lines of Fe i and 24 lines 
of Fe ii in the wavelength range 4400 - 8500 A. In spectra of the 
metal-poor stars most of these lines are very weak. Thus, 12 strong 
lines of Fe i were added to the analysis of the metal-poor stars. All 
line parameters are given in Tablesf2][3] and|4] 

The g/-values adopted in this work are weighted averages 
of different experimental valu es. An elaborate d iscussion of 
their accuracy can be found in IGehren et al. I d200ll) . The largest 



Oxford (Blackwell et al. 




19791; 


Blackwell, Petford & Shallis 


19791; Blackwell & Shallis 


1979 


Blackwell et al. 


198C; 


Blackwell, Shallis & SimmonsI 


198: 


; Blackwell et alj 


1982; 


Blackwell, Petford & SimmonsI 


198: 


; Blackwell et al. 1986). 



According to the NIST database, the uncertainties of the data 
are 3 - 10%. Smaller weights are assigned to th e g/- values of 
the H annover group dBard. Kock & Kock! Il99ll, iBard & Kockl 
1 1994 uncertainties typically 10 - 25%) and lO'Brian et alj dl99ll 
uncertainties 25 - 50%). The g/-values for the Fe n transitions 



7 Note that the uncertainty in mass is nominal and was determined by com- 
paring positions of the stars on the HRD with theoretical isochrones. 



Fe line formation for late-type stars 1 1 



Table 2. Parameters of the Fe i lines used for the solar analysis. 



Table 3. Parameters of the Fe n lines used for the solar analysis. 
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-31.4 


6978.85 


a 3 P 


z 3 P° 


1 


3 


2.48 


-2.480 


-31.4 


7401.68 


c 3 F 


w 3 D° 


5 


3 


4.19 


-1.500 


-31.2 


7912.87 


a 5 F 


z 7 D° 


11 


9 


0.86 


-4.848 


-31.9 


8293.51 


a 3 D 


y 3 D° 


5 


5 


3.30 


-2.203 


-31.3 



1 

,1 


Lower 


Upper 


Si 


go 


^low 


log?/ 


log Ce 


A 


level 


level 






|evj 






4491.40 


b 4 F 


z 4 F° 


4 


4 


2.856 


-2.71 


-32.1 


4508.29 


b A F 


z 4 D° 


4 


2 


2.856 


-2.44 


-32.0 


4576.34 


fe 4 F 


z 4 D° 


6 


6 


2.844 


-2.95 


-32.0 


4582.84 


b A F 


z 4 F° 


6 


8 


2.844 


-3.18 


-32.1 


4583.84 


b 4 F 


z 4 F° 


10 


8 


2.807 


-1.93 


-32.0 


4620.52 


b A F 


z 4 D° 


8 


8 


2.828 


-3.21 


-32.0 


4923.93 


a 6 S 


_6po 


6 


4 


2.891 


-1.26 


-32.1 


5018.44 


a 6 S 


z 6 P° 


6 


6 


2.891 


-1.10 


-32.1 


5169.03 


a 6 S 


z 6 P° 


6 


8 


2.891 


-1.00 


-32.1 


5197.58 


a A G 


Z 4 F° 


6 


4 


3.230 


-2.22 


-32.1 


5234.62 


a 4 G 


Z 4 F° 


8 


6 


3.221 


-2.18 


-32.1 


5264.81 


a 4 G 


z 4 D° 


6 


4 


3.230 


-3.13 


-32.0 


5284.11 


a 6 S 


_6po 


6 


8 


2.891 


-3.11 


-32.1 


5325.55 


a 4 G 


_4po 


8 


8 


3.221 


-3.16 


-32.1 


5414.07 


a 4 G 


Z 4 D° 


8 


8 


3.221 


-3.58 


-32.0 


5425.26 


a 4 G 


,4 F o 


10 


10 


3.199 


-3.22 


-32.1 


6239.95 


fo 4 D 


Z 4 P° 


2 


4 


3.889 


-3.41 


-32.0 


6247.56 


b 4 D 


.4po 


6 


4 


3.892 


-2.30 


-32.0 


6369.46 


a 6 S 


Z 6 D° 


6 


4 


2.891 


-4.11 


-32.1 


6432.68 


a 6 S 


z 6 D° 


6 


6 


2.891 


-3.57 


-32.1 


6456.38 


fo 4 D 


z 4 P° 


8 


6 


3.903 


-2.05 


-32.0 


6516.08 


a 6 S 


Z 6 D° 


6 


8 


2.891 


-3.31 


-32.1 


7222.39 


b 4 D 


Z 4 D° 


4 


2 


3.889 


-3.26 


-32.0 


7224.49 


fo 4 D 


z 4 D° 


2 


2 


3.889 


-3.20 


-32.0 


7515.83 


fo 4 D 


z 4 D° 


8 


6 


3.903 


-3.39 


-32.0 


7711.72 


b 4 D 


z 4 D° 


8 


8 


3.903 


-2.50 


-32.0 



Table 4. Parameters of the Fe I lines used for the abundance analysis of the 
metal-poor stars. g/-values are taken from NIST. 



,1 


Lower 


Upper 


Si 


Su 




log?/ 


logC 6 


A 


level 


level 






[eV] 






3581.19 


« 5 F 


z 5 G° 


11 


13 


0.86 


0.406 


-31.6 


3618.77 


a 5 F 


z 5 G° 


5 


7 


0.99 


-0.003 


-31.5 


3719.93 


a 5 D 


z 5 F° 


9 


11 


0.00 


-0.432 


-31.8 


3737.13 


fl 5 D 


z 5po 


7 


9 


0.05 


-0.574 


-31.6 


3745.56 


fl 5 D 


z 5po 


5 


7 


0.09 


-0.771 


-31.8 


3758.23 


a 5 F 


yV 


7 


7 


0.96 


-0.027 


-31.6 


3820.43 


a 5 F 


y 5 D° 


11 


9 


0.86 


0.119 


-31.6 


4045.81 


a 3 F 


y 3 F° 


9 


9 


1.49 


0.280 


-31.5 


4235.94 


; 7 D° 


e 7 D 


9 


9 


2.43 


-0.341 


-30.5 


4250.79 


fi 3 F 


z 3 G° 


7 


7 


1.56 


-0.714 


-31.5 


4415.12 


a 3 F 


z 5 G° 


5 


7 


1.61 


-0.615 


-31.5 


5586.76 


,5 F o 


e 5 D 


9 


7 


3.37 


-0.144 


-30.4 



are taken from iMelendez & Barbuvl j2009t). who renormalized 
the branching fractions from lRaassen & Uvlingsl ( fl998h and other 
theoretical s ources to the experimental lifetimes dSchnabel et al.l 
e.g., that of 1 19991 . |2004|) . To avoid biased results, we excluded 
any of their 'solar' values, which were obtained by a ID LTE 
spectroscopic analysis of the solar flux spectrum. That is, for the 
Fe ii transitions AA 5284.1,6239.95,6247.56,6456.38 A we kept 
the NIST-recommended values. 

Our adopted transition probabilites for the Fe i and Fe n lines 
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Figure 7. Comparison of transition probabilities adopted in this work (Ref) with that from the NIST database. 



are compared to the NIST-recommended values in Fig. [7] There 
are small differences for the weak Fe n transitions, which reflect 
the di screpancies between the N IST Fe n data based essentially 
on the | Raassen & Uvlingsl i ll 9981) stud y and the adopted va l ues of 
iMelendez & Barbuvl feOOa) . Note that iGrevesse & Sauvall fl999h 
in the re-analysis of the solar F e abundance discuss t he po ssibility 
of underestimated log gf 's fromE aassen & Uvlingsl < fl998h for the 
optical lines, although they do not find the same problem for the 
Fe n UV lines. 

Line broadening due to elastic collisions with H i atoms was 
com puted from the tables oflBarklem. Piskunov & O'Maral d2000h 
and iBarklem & Aspelund- Johansson! d2005t) . For the four Fe i 
lines (4445.47, 5600.22, 5661.35, 8293.51 A), broadening cross- 
sections a and their velocity dependence a were kindly provided 
by P. Barklem (2011, private communication). For seven lines, the 
lUnsoldl < ll955h values increased by a factor of 1.5 were adopted. In 
Tables|2][3] and and|4]these values are given in terms of commonly- 
used van der Waals damping constants log 



4.4 Sun and Procyon 

Owing to the superior quality of the solar spectrum, the initial spec- 
trum synthesis for the Sun was performed taking into account solar 
rotation with V rot G = 1.8 km s" 1 and a radial-tangential macrotur- 
bulence velocity £ RT , which was adjusted for each line separately. 
The typical values of £rt required to match the observed spectrum 
are 2.5 ... 4 km s _I for the profile fitting with the ID models, and 
1.6 ... 3 km s for that with t he (3D) model . We n ote that in full 
3D line formation calculations I Asplund et alj j2000h found that no 
macroturbulence was necessary to fit the line profiles due to the 
Doppler shifts from the convective motions, which are not explic- 
itly taken into account in the (3D) models used here. 

Selected Fe i and Fe n solar line profiles computed under 
LTE and NLTE with the (3D) models are compared with the ob- 
servations in Fig. [8] In general, NLTE Fe i line profiles are weaker 
compared to LTE, which is driven by the NLTE effect on the line 
opacity. For the strongest Fe i lines, forming very far out in the at- 
mospheres, deviation of the line source function from the Planck 
function, S ls < B V (T S ), leads to line core darkening. However, their 
wings, which dominate the total line strength, are formed at the 
depths where b t < 1 due to over-ionization, and the NLTE absorp- 
tion coefficient is smaller. Therefore the NLTE equivalent widths 



are smaller compared to LTE for a given abundance. As discussed 
in Sect. [3] NLTE effects on the Fe n lines are negligible. We also 
inspected the influence of inelastic H i collisions on the line profiles 
finding that the profile shapes are rather insensitive to the S H value. 
Any Fe line could be fitted by slightly adjusting the abundance or 
the damping parameters within the error bars (typically 10 percent). 

The mean Fe abundances computed with the mafags-os, 
marcs, and the (3D) model atmospheres are given in the Table [5] 
(see also Fig.ll2t, The total errors of the mean are only shown for 
the NLTE results obtained with the (3D) models. The total error 
was computed as: 



((c 



i)/2) 



1/2 



where, 



(Oie + °"tog£/,Fe + Ofogg^e + °"L.Fe)" 2 



The total error includes the following uncertainties: observational 
errors given by the standard deviation, errors stemming from the 
reference values of the surface gravity and effective temperatures 
(Table [T), and in the logg/ values. For the latter, we adopted the 
uncertainties given by NIST and assumed they are uncorrected, 
since the g/-values come from different sources in the literaturfl 
That is, the total uncertainty o"iog g / is decreased by a factor V^V, 
where N is the number of lines, and it takes into account the relative 
number of the Fe i and Fe n lines analyzed per star. 

We find that the excitation and ionization equilibrium of Fe 
in the Sun and Procyon is well satisfied under NLTE with both 
ID and (3D) models (Fig.[l2] Table[7]l. The statistical uncertainty 
of the mean abundance, which corresponds to the standard error 
of the line sample, is not larger than 0.01 dex. Some lines with 
W A > 100 mA seem to indicate slightly lower abundances, by 
~ 0.03 dex, compared to the rest of the sample both in ID and (3D) 
modelling. We do not assign much weight to these lines. First, it is 
almost impossible to discern the influence of abundance, damping, 
and weak blends in the line wings. Second, their asymmetric pro- 
files are clearly shaped by convective velocity fields, which are not 
accounted for in our models. As a consequence, it is hardly possi- 
ble to assign a unique 'best-fitting' abundance to such lines at all, 



Note that it is not possible to assess systematic errors in logg/'s. 
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Figure 8. Synthetic (black trace) and observed (dots) profiles of selected Fe i and Fe ii lines in the solar flux spectrum. The results are shown for the mafags-os 
solar model in NLTE (black trace) and LTE (dotted trace). 



that, in addition to a systematic error caused by the neglect of 3D 
effects, introduces subjectivity in abundance estimates. 

The mean NLTE abundances of Fe in the Sun and Procyon, 
determined using the (3D) model atmospheres, are 7.46 ± 0.02 re- 
spectively 7.43 ± 0.07 dex. ID models (marcs or mafags-os) yield 
slightly lower values, which are, however, all consistent within the 
total error. For the Sun, the latter only reflects the statistical uncer- 
tainty of the gf- values, as described above. We also find that adopt- 
ing LTE or scaling down the cross-sections for inelastic collisions 
with H i by an order of magnitude, Sh =0.1 does not yield any 
significant effect on the mean abundance. The difference between 
the two extreme cases A(NLTE-LTE) is rj 0.02 dex. Thus, based 



on the analysis of the solar-metallicity stars only it is not possible 
to single out the optimum value of Sh. It turns out, however, that 
this problem can be solved using the metal-poor stars with reliable 
stellar parameters (Sect. |4~5l . 

Comparing our results with the published values of the Fe 
abundance in the Sun and Procyon, we find full agreement with all 
data obtained with similar techniq ues and model atmospheres. We 
list only few most recent examples. iMelendez & Barbuvl fcOO^) us- 
ing a (3D) solar model very similar to that adopted in this work and 
L TE approach obtain lo g e Fe u = 7.45 dex. The NLTE ID estimate 
bv lGebjenetailfeOOll) obtained using the mafags-odf (note the dif- 
ference with the mafags-os version used here) models is log £pe,o = 
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7.48. ..7. 51 dex. In the follow-up NLTE study with the mafags-os 
models. IMashonkina et al.l d201ll ) find loge Fc = 7.41 . . . 7.56 de- 
pendi ng on the source of g f- va lues used for the Fe n lines. For Pro- 
cyon, IMashonkina et alj d20 1 lh claim discrepant abundances from 
the NLTE analysis of the Fe i and Fe n lines, log epci = 7.46 ± 0.07 
respectively log£ FcII = 7.52 ± 0.05 dex (their Table 3, solution 
Su = 0.1) suggesting that either increasing T cff by 80 K or de- 
creasing log g by 0.15 dex may solve the problem. Their results for 
Procyon with S H = 1, which is consistent with our study, are dis- 
crepant, log£p d = 7.41 ± 0.07 respectively log£ Fe n = 7.53 ± 0.05 
dex, which is most likely due to the different sources of log gf val- 
ues for the Fe n lines. 

We emphasize that further, more detailed, comparison to these 
and other studies is not meaningful because of various methodi- 
cal differences, such as atomic data, line selection, and microturbu- 
lence, to name just a few. 



4.5 Analysis of metal-poor stars 

We start with a discussion of how NLTE affects the Fe abundances 
obtained with ID and (3D) models, delineating the key differences. 
In Sect. 14.5.21 we then illustrate how these differences propagate in 
the determination of effective temperatures and surface gravities. 



4.5.1 Metallicities 

The most important test of our new models is whether they are able 
to recover ionization equilibrium of Fe i and Fe n for the stars with 
parameters determined by independent techniques, given a unique 
Sh. As described in Sect. 14. 21 a small set of metal-poor stars in dif- 
ferent evolutionary stages was selected for this purpose with T e $ 
known from the interferometry and/or IRFM, and gravities well 
constrained by parallaxes. Metallicities and microturbulence pa- 
rameters were then determined keeping the T c jf and log g fixed, and 
varying the value of Sh in the NLTE calculations. 

We find that the solution providing optimum ionization bal- 
ance of Fe for all metal-poor stars is achieved with Sh = 1, i.e, un- 
sealed Drawin's H i collisions (see also discussion in Sect. 14.5.21 . 
The Fe abundances determined in this way, along with the opti- 
mized (3D) NLTE microturbulence £ t , are given in the Table [5] 
These are mean abundances from the Fe i and Fe n lines, and 
the errors were computed as described in Sect. 14.41 For all stars, 
but the Sun and G 64-12 with its very uncertain parallax, the er- 
rors are dominated by the uncertainties of the reference T^g. The 
abundances are shown in Figs.[l2][l3] and!14las a function of line 
equivalent widths and lower level excitation potentials (in eV). The 
error bars in the figures indicate the size of NLTE abundance cor- 
rections for the low-excitation (< 2.5 eV) and high-excitation Fe i 
lines. They are levelled out at the mean NLTE abundance obtained 
from the Fe i and Fe n lines with Sh = 1. The bar's upper and 
lower ends correspond to the mean NLTE abundance obtained with 
Sh =0.1 and the LTE abundance, respectively. The abundances av- 
eraged over the measured Fe i respectively Fe n lines are presented 
in Table [7] along with their one sigma errors. Note that these are 
not total propagated errors as described in Sec. 4.4, but standard 
deviations of the mean abundance for each ionization stage. These 
quantities mainly serve to show the difference in the line scatter 
between the ID (marcs) and (3D) models. 

The major differences between ID and (3D), as well as LTE 
and NLTE, are apparent from inspection of the Fe abundances in the 



Tableland the NLTE abundance correction^] shown in Fig. [9] For 
the Sun and Procyon, the NLTE effects are rather subtle. The differ- 
ence between the LTE and NLTE abundances determined from the 
Fe i lines is even slightly smaller in the (3D) case compared to ID 
(Fig-IU top panels). The basic reason is that the temperature gra- 
dient in the solar-metallicity marcs models is slightly steeper than 
in the (3D) models (Fig. [TJ. At logT 50 oo > 0, where the UV-blue 
continuum is formed, (dT/dT) lD > (dT/dr) 3D , whereas in the outer 
layers the gradients are very similar. Over-ionization in Fe i by ra- 
diation field emerging from the hot deep layers is, thus, somewhat 
stronger in the marcs models. 

Over-ionization is more important at low metallicities (see 
Sect. [3] of this paper, and a more extensive discussion in Paper II). 
This causes non-negligible differences between the NLTE and LTE 
abundances inferred from the Fe i lines. The differences are in the 
range 0.05 - 0.15 dex for ID models and weak Fe i lines, but in- 
crease to 0.2 - 0.5 dex for the strong low-excitation Fe i lines with 
the (3D) models (Fig. [9}. The larger NLTE effects of metal-poor 
(3D) models compared to ID models are due to the larger decou- 
pling between the radiation temperature and the local temperature, 
owing to the dramatically cooler outer layers of the (3D) models 
compared to ID (see Sect.[3]and Fig.[6](. NLTE corrections to the 
Fe ii lines are negative and are not larger than -0.03 dex. As a re- 
sult, for the metal-poor stars, the difference between the mean LTE 
and NLTE Fe abundances is, at least, twice as large in the (3D) 
case compared to that obtained with marcs and mafags-os models. 
On the other side, the same phenomenon, which leads to ampli- 
fied NLTE effects in the (3D) calculations, i.e., very steep (dT/dr) 
above logrsooo ~ 0, is also responsible for the LTE strengthen- 
ing of the Fe i lines in (3D) compared to ID. The Fe n lines are 
only slightly weaker in (3D) calculations. Our LTE ID metallicities 
are, thus, systematicall y larger than the (3D) results, in agree ment 
with other studies (e.g., ICollet. Asplund & Tramp edach 2007]). We 
note , however, that dire ct comparison of our results with that of, 
e.g.. ICollet et al.l J2007I) . is not meaningful since they performed 
LTE line formation calculations with full 3D hydrodynamic model 
atmospheres. 

The fact that the (3D) models also predict somewhat lower 
LTE Fe abundances, yet larger NLTE abundance corrections than 
ID models, has the important consequence that the mean NLTE 
metallicities in ID and (3D) turn out to be in agreement for the 
whole range of stellar parameters (Alog£p c < 0.1). 

However, even though the mean metallicities are similar, indi- 
vidual Fe i line abundances are still somewhat discrepant, and show 
systematic trends with line equivalent width and excitation poten- 
tial of the lower level of a transition, E\ ow (Figs. [T3] I14t . NLTE 
abundances inferred from the low-excitation lines of Fe i with ID 
models are too large compared to the high-excitation lines and lines 
of Fe n. A very similar picture is obtained with the LTE approxi- 
mation, and it reverses sign in LTE (3D). A combination of (3D) 
and NLTE alleviates this discrepancy, so that Fe i and Fe n abun- 
dances become more consistent. The improvement comes from the 
different sensitivity of the NLTE effects to E low in the ID and (3D) 
case. This is clearly seen in Fig. [9] A least-square fit to the NLTE 
abundance corrections is also shown, which is a simple albeit very 
illustrative measure of the mean NLTE effect on the excitation bal- 
ance of Fe i. In the ID case, the Fe i lines have rather similar NLTE 



9 NLTE abundance correction, Ap e = log A(Fe) NLTE - logA(Fe) LTE , is the 
difference between the NLTE and LTE abundances required to fit a line with 
a given equivalent width 
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Table 5. Metallicities [Fe/H] and microturbulence £ t (km s ) determined relative to the Sun for the reference stars. The g t for the mafags models are identical 
to that obtained with the marcs models and are not shown. NLTE results refer to 5h = 1 . The references to the T^g and log g values are given in the Table 1 . 
See text. 



Star fixed derived — » 

T c ff logg MARCS MAFAGS (3D) ITtot 

NLTE £ LTE £ NLTE LTE NLTE £ LTE f NLTE, (3D) 



Sun 


5777 


4.44 


7.44 


1.02 


7.43 


1.02 


7.46 


7.45 


7.46 


1.08 


7.45 


1.09 


0.02 


Procyon 


6543 


3.98 


-0.06 


1.93 


-0.08 


1.94 


-0.05 


-0.07 


-0.03 


2.05 


-0.04 


2.06 


0.07 


HD 84937 


6408 


4.13 


-2.07 


1.42 


-2.10 


1.43 


-2.05 


-2.09 


-2.03 


1.38 


-2.12 


1.35 


0.08 


HD 140283 


5777 


3.70 


-2.38 


1.26 


-2.42 


1.27 


-2.40 


-2.43 


-2.40 


1.18 


-2.53 


1.17 


0.07 


HD 122563 


4665 


1.64 


-2.55 


1.60 


-2.60 


1.61 


-2.48 


-2.52 


-2.57 


1.66 


-2.66 


1.67 


0.24 


G 64-12 


6464 


4.30 


-3.12 


1.75 


-3.24 


1.75 


-3.16 


-3.23 


-3.10 


1.75 


-3.37 


1.75 


0.17 



Table 6. Spectroscopic 7" c ff, logg, g t (km s ), and [Fe/H] values. The parameters, which were kept fixed in the calculations, are given in the column 2. 



Star fixed derived — » 

log^eff ID, LTE ID, S H = 1 ID, S H = 0.1 (3D), LTE (3D), S H = 1 (3D), 5 H = 0.1 



Sun 


4.44 


5794 


1.0 


7.44 


5780 


1.0 


7.45 


5753 


1.0 


7.46 


5805 


1.1 


7.47 


5793 ± 35 


1.1 


7.47 ± 0.03 


5775 


1.0 


7.48 


Procyon 


3.98 


6631 


1.9 


-0.04 


6578 


1.9 


-0.05 


6487 


1.9 


-0.05 


6576 


2.0 


-0.04 


6546 ± 43 


2.0 


-0.04 ± 0.04 


6499 


2.0 


-0.05 


HD84937 


4.13 


6405 


1.4 


-2.11 


6328 


1.4 


-2.13 


6115 


1.4 


-2.16 


6483 


1.3 


-2.08 


6314 ±69 


1.4 


-2.10 ± 0.07 


6134 


1.5 


-2.12 


HD 140283 


3.70 


5772 


1.2 


-2.43 


5703 


1.2 


-2.45 


5510 


1.2 


-2.49 


6078 


1.2 


-2.36 


5820 ± 65 


1.2 


-2.38 ± 0.07 


5586 


1.3 


-2.41 


HD122563 


1.64 


4754 


1.6 


-2.51 


4705 


1.6 


-2.51 


4609 


1.5 


-2.49 


4838 


1.7 


-2.53 


4755 ± 59 


1.7 


-2.50 ± 0.14 


4669 


1.7 


-2.48 


G64-12 


4.30 


6592 


1.8 


-3.15 


6471 


l.S 


-3.17 


6006 


1.8 


-3.34 


6757 


1.8 


-3.18 


6439 ± 152 


1.8 


-3.17 ±0.17 


6258 


1.8 


-3.16 


Sun 


5777 


4.40 


1.0 


7.44 


4.43 


1.0 


7.45 


4.48 


1.0 


7.47 


4.38 


1.0 


7.46 


4.40 ± 0.06 


1.0 


7.47 ± 0.03 


4.44 


1.0 


7.48 


Procyon 


6543 


3.81 


1.9 


-0.11 


3.93 


1.9 


-0.07 


4.08 


1.9 


-0.04 


3.92 


2.0 


-0.05 


3.99 ± 0.14 


2.0 


-0.04 ± 0.08 


4.06 


2.0 


-0.02 


HD84937 


6408 


4.14 


1.4 


-2.11 


4.26 


1.4 


-2.07 


4.57 


1.4 


-1.98 


4.00 


1.3 


-2.12 


4.28 ± 0.16 


1.4 


-2.04 ± 0.08 


4.53 


1.5 


-1.95 


HD 140283 


5777 


3.72 


1.2 


-2.42 


3.84 


1.2 


-2.39 


4.21 


1.2 


-2.26 


3.11 


1.1 


-2.59 


3.63 ± 0.14 


1.1 


-2.41 ± 0.08 


4.05 


1.3 


-2.26 


HD 122563 


4665 


1.22 


1.6 


-2.65 


1.48 


1.6 


-2.57 


1.83 


1.6 


-2.45 


0.91 


1.7 


-2.74 


1.31 ± 0.26 


1.7 


-2.59 ± 0.19 


1.61 


1.7 


-2.48 


G64-12 


6464 


4.06 


1.8 


-3.26 


4.29 


1.8 


-3.18 


4.90 


1.8 


-3.05 


3.75 


1.8 


-3.37 


4.34 ± 0.20 


1.8 


-3.16 ± 0.10 


4.56 


1.8 


-3.04 



corrections irrespective of their excitation potential and equivalent 
width. In the (3D) case, however, the low-excitation transitions, 
E\ ovl < 2 eV, tend to experience significantly larger departures from 
LTE compared to the higher-excitation transitions. The reason is 
that, owing to the dramatically cooler surfaces of the metal-poor 
(3D) models, low-excitation Fe i lines become stronger and more 
sensitive to the physical conditions in the outer atmospheric lay- 
ers, where the influence of non-local radiation field is extreme, not 
only decreasing the line opacity but also pushing the line source 
functions to super-thermal values. Thus, the NLTE (3D) profiles of 
these lines are weaker compared to LTE. The magnitude of this ef- 
fect depends on stellar parameters (see also the discussion in Sect. 
[3j. As a consequence, one would also expect that excitation bal- 
ance achieved under LTE with the (3D) models will strongly over- 
estimate T e g, and the error shall increase with T e g and decreasing 
metallicity as indicated by the increasing slope of A(NLTE - LTE) 
vs £■[„„ for more metal-poor and hotter stars. 

A corollary is that accurate metallicities for late-type stars can 
be obtained with standard ID model atmospheres if NLTE effects 
in Fe i are taken into account. However, this is true only if the fol- 
lowing condition is satisfied: a sufficiently large number of Fe i 
and Fe n lines of different types are included in an analysis, so that 
individual line-to-line abundance discrepancies cancel out. This ap- 
proach, however, does not eliminate residual trends of abundance 
with line equivalent width and excitation potential. Although these 
trends can be, at least in part, corrected for by adjusting microtur- 



bulence, a better approach is to restrict an analysis with ID hydro- 
static models to hig h-excitation Fe i l ines only, as also advocated in 
the literature (e.g., iGustafssonll 19831) . 

4.5.2 Effective temperatures and surface gravities 

The spectroscopic effective temperatures and surface gravities de- 
termined according to the procedure described in Sect. |4. H are given 
in the Table|6] 

Fig. [10] (top panel) plots the effective temperatures obtained 
for the reference stars adopting a fixed surface gravity inferred from 
stellar parallaxes. The results are shown for the three cases: LTE, 
NLTE with the standard model atom (Sh = 1), and NLTE with 
Drawin's H i collision cross-sections scaled by a factor of 5h = 0. 1. 
The reference temperature scale, T eff REF , refers to the IRFM results 
(see Sect, \4.2l for all stars but Procyon with the interferometric 
angular diameters; their uncertainties are shown with blue filled ar- 
eas. The uncertainties of our spectroscopic r cff 's were estimated 
by mapping the errors in atomic transition probabilities, observa- 
tional errors, and errors in the input surface gravity into a range in 
temperatures. In order not to overload the figures, these errors are 
shown for the NLTE 5h = 1 solution only. We emphasize that the 
7Ws were obtained by enforcing ionization balance, but not requir- 
ing null slopes with excitation potential for the reasons discussed 
above. 

The overall best result in terms of consistency between the 
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Table 7. LTE and NLTE abundances determined from the Fe I and Fe II lines using the reference values of r c ff and logg (Tables 1 and 5). NLTE results refer 
to Sh = 1 ■ Standard deviations cr of the Fe I respectively Fe II mean abundances are also indicated. 

Star marcs mafags (3D) 

LTE NLTE LTE NLTE LTE NLTE 

Fel Fell Fe I <r Fell cr Fe I Fell Fe I Fell Fe I Fell Fe I cr Fell cr 



Sun 


7.43 


7.45 


7.44 


0.05 


7.44 


0.04 


7.46 


7.44 


7.48 


7.43 


7.44 


7.47 


7.45 


0.04 


7.47 


0.04 


Procyon 


7.33 


7.40 


7.37 


0.04 


7.40 


0.04 


7.39 


7.38 


7.43 


7.37 


7.40 


7.43 


7.43 


0.04 


7.43 


0.04 


HD 84937 


5.32 


5.33 


5.39 


0.07 


5.33 


0.04 


5.36 


5.38 


5.42 


5.38 


5.30 


5.38 


5.45 


0.06 


5.38 


0.03 


HD 140283 


5.02 


5.01 


5.07 


0.09 


5.01 


0.04 


5.03 


5.01 


5.09 


5.01 


4.86 


5.09 


5.05 


0.05 


5.09 


0.04 


HD 122563 


4.78 


4.95 


4.87 


0.12 


4.95 


0.05 


4.93 


4.95 


5.00 


4.95 


4.71 


5.01 


4.85 


0.08 


5.01 


0.05 


G 64-12 


4.25 


4.23 


4.36 


0.05 


4.24 


0.09 


4.23 


4.21 


4.34 


4.22 


4.20 


4.28 


4.39 


0.05 


4.30 


0.08 
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Figure 9. The NLTE abundance corrections to the Fe i lines obtained with Marcs and (3D) model atmospheres for the same reference T e g and log g (see text). 
The metallicities indicated are that obtained from the NLTE spectroscopic analysis with the (3D) model atmospheres. 
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Figure 10. Comparison of spectroscopic stellar parameters, r c ff and logg, with the data obtained by independent methods, IRFM and parallaxes. G 64-12 is 
not shown in the bottom figure, because its astrometric gravity is too uncertain. 



spectroscopic effective temperatures, T c ff_Fc, and effective temper- 
atures obtained by independent methods, T e ff_REF, is achieved un- 
der NLTE using the reference model atom with Sh = 1, i.e. stan- 
dard Drawin's H i collisions. The ID and (3D) models perform 
similarly, predicting only minor offsets of T c ffFe from the refer- 
ence IRFM r e (f scale, with an rms error of the order of ~ 50 K. 
Assuming LTE or decreasing the efficiency of inelastic H i colli- 
sions in the NLTE statistical equilibrium calculations by an order 
of magnitude, Sh = 0.1, leads to larger discrepancies. The latter 
result is very important. At low metallicity NLTE effects in Fe i 
are extremely sensitive to the magnitude of collisional thermaliza- 
tion by H i. Our analysis, thus, not only confirms that it is possi- 
ble to constrain the absolute magnitude of atomic collision cross- 
sections using metal-poor stars, but also demonstrates that a single 
scaling factor (unity in our case) to the classical Drawin's formula 
is a reasonable approximation. What concerns the LTE T eff values, 
they are still within the errors of T cff REF in the ID case. We cau- 
tion, however, that the NLTE effects for our stars are small in ID, 
and this statement may not apply for stellar parameters where the 
NLTE effects are more significant, i.e. for stars with higher T cB and 
lower logg. The discrepancy between the LTE spectroscopic and 
the IRFM scale increases to 200 K for the metal-poor stars when 
(3D) models are employed, which is consistent with the discussion 
in the previous section. 

Our second approach, in which the r cl i values are kept fixed 
and surface gravities are obtained by forcing Fe ionization equi- 
libria, returns similar results. Fig. [TO] (bottom panel) compares the 
spectroscopic logg's with the values determined from parallaxes, 
the shaded regions indicate the uncertainty of the latter (see Sect. 



\4.2\ and the error bars of the former account for the uncertainties 
in logg/, IRFM effective temperatures, and observational uncer- 
tainties. As seen from this figure, only the NLTE approach with 
Sh = 1 gives log g values on average consistent with that inferred 
from parallaxes. In (3D), LTE gravities are systematically too low. 
In ID, the discrepancy is smaller, but appears to increase with T e g. 
The assumption of Sh = 0. 1 brings the Fe i abundances so high, 
that much larger gravities are generally needed to match Fe i with 
Fe n. 



4.5.3 Comparison with other studies 

Korn et al. (2003) performed a study similar to ours aiming to con- 
strain the efficiency of inelastic H i collisions from the NLTE ID Fe 
calculations for a small sample of late-type stars. Using the mafags- 
os model atmospheres and the Drawin's formulae for H i collision 
cross-sections, they determined Sh = 3 as a scaling factor to the 
latter yielding satisfactory ionization balance of Fe i and Fe n. For 
the stars in common, HD 140283, HD 84937, and Procyon, their 
results are consistent with our values for Sh = 1. This study was, 
however, later superseded b y that of Mashonkina et al. (2011), in 
which the Fe model atom o f iGehren et ail feOOlL b), also used by 
Korn et al. (2003), was updated by more recent atomic data. This 
work, in fact, showed that a lower efficiency of inelastic H i col- 
lisions (Sh = 0.1) is needed to satisfy ionization balance of Fe. 
Their results for HD 84937 and HD 122563 are in agreement with 
our values within the error bars: -2.00 ± 0.07 (Fe I) respectively 
-2.08 ± 0.04 (Fe II) for HD 84937, and -2.61 ± 0.09 (Fe I) re- 
spectively -2.56 ± 0.07 (Fe II) for HD 122563. They also find that 
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the NLTE Fe i-based abundances are correlated with line excita- 
tion potential, with a slope of 0.013 dex eV~' for HD 84937 and 
-0.054 dex eV~' for HD 122563. From Fig. [T3] we can not con- 
firm their slope for HD 84937 with ID models, however, we would 
get exactly same results, adopting Mashonkina et al. linelists. Their 
positive slope for HD 84937 would likely had vanished had they 
included few more subordinate near-UV Fe i lines with very accu- 
rate gf- values in their analysis. For HD 122563, however, we also 
find a small discrepancy between the high-excited Fe i and Fe n 
lines, in addition to residual slopes of the NLTE Fe abundances 
with E\ ow : -0.11 dex eV~' for the mafags-os, -0.09 dex eV~' for 
marcs, and -0.05 dex eV~' for the (3D) model. One of the impor- 
tant differences between our and their study, which also explains 
our more extreme slopes in ID, is the way we handle microturbu- 
lence. We use Fe n lines, whereas Mashonkina et al. (2011) rely on 
the Fe i lines. Nevertheless, the fact that our results obtained with 
ID model atmospheres support the findings of Mashonkina et al. 
(201 1) is reassuring. 



4.5.4 HD 122563 

HD 122563 is one of the most well-studied metal-poor giants. Nev- 
ertheless, its parameters are badly constrained. For example, Ful- 
bright (2000) gives T cff = 4425 K, logg = 0.6, and [Fe/H] = -2.6, 
whereas according to Mashonkina et al. (201 1) the star is better de- 
scribed by a model with r cff = 4600 K, logg = 1.6, and [Fe/H] 
~ -2.58. 

At Sh = 1, no combination of parameters can achieve a fully 
satisfactory spectroscopic solution for this star, consistently achiev- 
ing adequate excitation and ionization balance. This is evident from 
the residual offset between the high-excitation Fe i and Fe n lines 
in Fig. [14] In NLTE with Sh = 1, the spectroscopic logg is 0.3 
dex lower than the value derived from the Hipparcos parallax. The 
discrepancy can be minimized by adopting less efficient hydrogen 
collisions, which, however, contradicts the spectroscopic results of 
the other metal-poor stars. 

Based on HD 122563 alone, the evidence for a lower Sh is thus 
limited. Similar investigation of larger samples of metal-poor giants 
and dwarfs with well-defined parameters may help to better clarify 
the situation. 



4.6 Comparison with evolutionary tracks 

In an attempt to assess the trustworthiness of the spectroscopic re- 
sults versus those obtained from independent measurements, we 
compared them to evolutionary tracks. The position of the four 
metal-poor stars in the r cff -logg plane is shown in Fig. [TT] The 
black cross with error bars indicated is the result with the IRFM 
T e g and astrometric surface gravity, and the blue square and red di- 
amond correspond to our spectroscopic (3D) NLTE r eff and logg, 
respectively. For comparison, evolutionar y tracks with close pa- 
rameters computed with the garstec code dWeiss & Schlattl|[2008h 
are shown. These were kindly made available to us by A. Serenelli 
(201 1, private communication). In order not to overload the figures, 
we show only the NLTE (3D) spectroscopic results, while briefly 
discussing other solutions from the Table [6] below. 

For all metal-poor stars, a combination of astrometric gravity 
and spectroscopic r c ff/[Fe/H], obtained with NLTE Sh = 1 or LTE 
with ID and (3D) models, appears to be most consistent with stel- 
lar evolution calculations. The error bars of our data are omitted in 
the plots for clarity (they are given in the Table |6j. ID LTE r^'s 



are somewhat larger than the ID NLTE results, that makes the stars 
a bit more massive. It is interesting that even with (3D) LTE ef- 
fective temperatures, which are ~ 50 - 300 K larger compared to 
NLTE results, the stars still match the tracks with adequate mass 
and metallicity, which, in fact, also implies more realistic ages. 

A combination of IRFM T eff and the spectroscopic (3D) LTE 
log g/[Fe/H] makes the agreement with stellar evolution less satis- 
factory. In particular, HD 140283 is now in the 'forbidden' zone of 
the HRD. Situation is even worse for the NLTE Sh = 0. 1 results, 
obtained with ID and with (3D) model atmospheres. In this case, 
all four metal-poor stars are some ~ 2 - 3cr in and log g away 
from any reasonable track. 

HD 122563 is a somewhat outstanding object. Assuming the 
tracks to be realistic, it also becomes clear that the IRFM T s g and/or 
luminosity of the star inferred from its relatively accurate parallax 
are too low. Other attempts, such as variation of the mass or de- 
creasing the ^-enhancement, either do not have any effect or make 
the discrepancy with masses and ages even larger. In this respect, 
it is encouraging that our spectroscopic (3D) NLTE values of T eff 
and logg bring the star into a much better agreement with stellar 
evolution, although these value appear to be shifted away from the 
presumably more reliable IRFM temperature and astrometric grav- 
ity. 



5 SUMMARY AND CONCLUSIONS 

One important problem in spectroscopy of late-type stars is to quan- 
tify systematic errrors in the determination of Tgg , log g, and [Fe/H] 
from the Fe i/Fe n excitation-ionization equilibrium due to the as- 
sumptions of LTE and hydrostatic equilibrium in the models. 

We have constructed a new extensive NLTE Fe model atom 
using the most up-to-date theoretical and experimental atomic data 
available so far. The model has been applied to the analysis of the 
Sun and a number of well-studied late-type stars with parameters 
determined by other independent methods. We have used standard 
ID LTE marcs and mafags-os model atmospheres, as well as aver- 
age stratifications computed from full time-dependent, 3D, hydro- 
dynamical simulations of stellar surface convection. In addition, 
we compared the commonly-used NLTE line formation program 
packages detail/siu and multi. Despite obvious differences in their 
numerical scheme and input physical data, we find that the final re- 
sults are consistent both in terms of the NLTE statistical equilibria 
and absolute LTE and NLTE abundances. 

Our LTE and NLTE results for the ID models are similar to 
most of the previous findings in the literature. Statistical equilib- 
rium of Fe, which is a minority ion at the typical conditions of 
these cool and dense atmospheres, favors lower number densities 
of Fe i compared to LTE. The number densities of Fe n are hardly 
affected by NLTE. In general, this leads to a weakening of Fe i lines 
compared to LTE, which, in turn, requires larger Fe abundance to 
fit a given observed spectral line. The magnitude of departures from 
LTE depends on stellar parameters. With ID hydrostatic model at- 
mospheres the NLTE corrections on Fe i lines do not exceed 0. 1 dex 
for stars with [Fe/H] > -3.5, and they are negligible for the Fe n 
lines. The situation changes dramatically for the mean 3D ((3D)) 
models with their cooler surfaces, and, thus, steeper T(r) relations. 
The NLTE abundance corrections can be as large as ~ 0.5 dex for 
the resonance Fe i lines formed in very metal-poor atmospheres. In 
contrast to LTE, NLTE strengths of the Fe i lines predicted by the 
(3D) and ID models are rather similar, because the line formation 
is largely dictated by the radiation field forming around the opti- 
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Figure 11. The position of the four metal-poor stars in the Teff -log g plane. The black cross with error bars is the result with the IRFM 7" e ff and astrometric sur- 
face gravity, and the blue square and red diamond correspond to the spectroscopic NLTE (3D) T s g and log g from this work, respectively, garstec evolutionary 
tracks are also shown for comparison. For G 64-12 only spectroscopic logg/IRFM r c ff solution is shown. 



cal surface, log tjooo ~ 0, where the thermal structures of (3D) and 
ID models are similar. This suggests that the full 3D NLTE results, 
once they become feasible, will be even closer to our (3D) NLTE 
calculations. 

The solar Fe abundance obtained under NLTE using the (3D) 
model atmospheres is 7.46 ± 0.06 dex, and the uncertainty is fully 
dominated by the errors of the experimental transition probabilities. 
Adopting LTE or ID hydrostatic model atmospheres lowers the so- 
lar Fe abundance by ~ 0.02 dex. For Procyon we infer slightly 
sub-solar metallicity, [Fe/H] = -0.03 dex. 

We find that the Fe r/Fe n ionization balance can be well es- 
tablished for all reference stars with the ID and (3D) model at- 
mospheres and the NLTE model atom with Drawin's H i collision 
cross-sections. Strong resonance and subordinate Fe i lines are very 
sensitive to the atmospheric structure, thus, classical ID models fail 
to provide consistent excitation balance. A better agreement among 
Fe i lines spanning a range of excitation potential is obtained with 
the mean 3D models, although the optimum solution, which fully 
eliminates the correlation of the abundance and line excitation po- 
tential, may necessitate self-consistent NLTE modelling with full 
3D hydrodynamical model atmospheres. The assumption of LTE in 
combination with (3D) models leads to large errors in T e g and log g 
inferred from the Fe ionization balance, yet we find satisfactory re- 
sults for certain combinations of stellar parameters using standard 
hydrostatic model atmospheres. 

The mean NLTE metallicities determined with ID and (3D) 



models are in agreement in the range of stellar parameters investi- 
gated here, although there are marked residual slopes of abundance 
with line excitation potential and strength. We can, thus, conclude 
that accurate, albeit not precise, metallicities for late-type stars can 
be obtained with classical hydrostatic ID model atmospheres if 
NLTE effects in Fe i are taken into account and a sufficiently large 
number of Fe i and Fe n lines with different excitation potentials is 
used, so that individual line-to-line abundance discrepancies can- 
cel out. An alternative solution, which avoids additional biases in- 
troduced by adjusting microturbulence parameter, is to use high- 
excitation Fe i lines only. Thus, the results obtained with classical 
ID models depend on the choice of the line list and, in particular, 
on the balance of the number of lines of different types. A combi- 
nation of (3D) and NLTE alleviates this problem. 

Our NLTE effective temperatures and gravities are consis- 
tent with the parameters determined by other less model-dependent 
methods, in particular with IRFM r c ff 's and logg's inferred from 
parallaxes. The results for the metal- poor dwarfs and subgiant are 
also consistent with stellar evolution predictions. For HD 122563, 
the results are inconclusive. In addition to a residual slope of Fe i- 
based abundances with line excitation potential, the NLTE (3D) 
spectroscopic gravity is roughly 0.3 dex lower than the astromet- 
ric result. On the one side, the latter is somewhat uncertain due to 
the difficulty of determining the mass of this metal-poor giant. Fur- 
thermore, comparisons with evolutionary tracks favor higher T tf! 
and lower logg which would be more consistent with our spec- 
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troscopic estimates. On the other side, the discrepancies might be 
indicative of remaining systematic uncertainties in the stellar atmo- 
sphere models, such as neglect of convective inhomogeneities and 
atmospheric extension, which become increasingly important with 
increasing and decreasing log g. The importance of the latter is 
also amplified at higher metallicities. 
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Table 1. Equivalent widths of the Fe i lines (in units of mA). 



Line (A) 
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Procyon 


HD 84937 


HD 140283 


HD 122563 
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3737.13 
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3745.56 
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3758.23 












64.00 


3820.43 
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70.40 


4045.81 






115.95 






71.02 


4235.94 






52.55 


59.23 


102.58 


10.30 


4250.79 
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78.99 


128.14 


24.26 


4415.12 






75.30 


84.90 
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8.49 
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34.62 


83.36 
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111.49 


87.12 


7.33 


13.49 


79.57 




5044.21 


77.48 
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2.15 


9.35 
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101.08 


75.32 


4.25 
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129.41 
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20.41 


79.26 




5225.53 


73.14 
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Figure 12. NLTE abundances determined from the low-excitation Fe i (filled circles), high-excited Fe i (blue diamonds), and Fe n (squares) lines using 
marcs-os, mafags-os, and (3D) model atmospheres for the Sun (top) and Procyon (bottom). Small vertical bars on the left and right-hand side of each plot 
show the size of NLTE abundance corrections for the low-excited (< 2.5 eV) and high-excited Fe i lines, respectively. They are leveled out at the mean NLTE 
abundance obtained from the Fe i and Fe n lines with 5h = 1 ■ The bar's upper and lower ends correspond to the mean NLTE abundance obtained with Sh = 0. 1 
and the LTE abundance, respectively. The dark and light shaded regions indicate the standard and total errors of the mean (see text). 
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Figure 13. Abundances determined from the Fe I and Fe II lines using stellar parameters from Table[5]for HD 84937 (top) and HD 140283 (bottom). 
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Figure 14. Abundances determined from the Fe I and Fe II lines using stellar parameters from Table[5]for HD 122563 (top) and G 64-12 (bottom). Note that 
the total errors are large as they also reflect the uncertainties in the input stellar parameters. 
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Table 2. Equivalent widths of the Fe n lines (in units of mA). 



Line (A) Sun Procyon HD 84937 HD 140283 HD 122563 G 64-12 
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